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ABSTRACT 

A photon of wavelength A ~ 1 [Am interacting with a dust grain of radius a p ~ 1mm (a "pebble") 
undergoes scattering in the forward direction, largely within a small characteristic diffraction angle 
6 S ~ X/a p ~ 100". Though mm-size dust grains contribute negligibly to the interstellar medium's 
visual extinction, the signal they produce in scattered light may be detectable, especially for variable 
sources. Observations of light scattered at small angles allows for the direct measurement of the large 
grain population; variable sources can also yield tomographic information of the interstellar medium's 
mass distribution. The ability to detect brilliant pebble halos require a careful understanding of the 
instrument PSF. 
Subject headings: dust scattering 
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1. INTRODUCTION 

The size distribution of dust in the diffuse interstellar 
medium has traditionally been studied through the wave- 
length dependence of interstellar extinction and polariza- 
tion, and o bservations o f scattered light from reflection 
nebulae (cf. iDraine 200 3]). These studies have resulted in 
models (fWeingartner fc Draind l200H: IZubko et al1l2004l : 
IDraine fc Fraissdl2009f) ." where the bulk of the dust mass 
resides in grains with radii a < 0.5 /xm, with approxi- 
mately 50% of the grain mass above and below a char- 
acteristic radius ~ 0.15 i/m. To account for the observed 
extinction, these models tend to consume the bulk of 
available elements such as Mg, Si, and Fe in the "ob- 
served" grain population with a < 0.5 /Ltm. Little is 
known regarding the dust grain size distribution above 
a grain size of order ~ 0.5 ^m, aside from the general 
expectation that grains above this size should contain at 
most a small fraction of the total dust mass. 

Therefore, it was surprising when impact detectors 
on Ulysses and Galileo measure d a flux of dust par- 
ticles entering th e heliosphere (jLandgraf et al.l 120001 : 
iKruger et a,l.ll2007h that appeared to indicate that the in- 
terstellar grain size distribution extended to much larger 
grains, with approximately equal mass per unit logarith- 
mic interval out to the largest sizes (awl /im) that 
could be detected. This finding was completely at odds 
with the co nclusion s draw n from studies of interstellar 
extinction. iDraind (|2009h argued that the size distri- 
bution inferred for the dust particles approaching the 
heliosphere could not characterize average interstellar 
dust, but the situation remains unclear. It is further 
confounded by radar detection of a ~ 30 /xm particles 
entering the Earth's atmosphere on solar-hyperb olic tra- 
jectories (jTavlor et al.lll99alBaggalgy||2000l . l200l imply- 
ing that they are arriving from interstellar space. The 
mass flux in these particles exceeds the mass flux in 
0.2 < a < 1 /xm particles inferred from the Ulysses and 
Galileo observations. If these particles are truly entering 
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the solar system from interstellar space, it implies a total 
mass density in a > 0.5/im grains that, at least locally, 
exceeds that in the a < 0.5 /xm grains in conventional 
models for interstellar grains. Such an abundance of very 
large grains would be difficult to understand given the 
limitations of interstellar abundances of elements from 
which such grains would be constituted. 

In this Letter, we show that very large interstellar 
grains, or "pebbles," are detectable in scattered light, 
a phenomenon we refer to as "brilliant pebbles." The 
basic idea is presented in |j2j Then, in <j3]we assess the 
observability of these scattered light halos, and show that 
detection is possible even if only a few percent of the in- 
terstellar grain mass is in mm-sized particles. Variable 
sources (e.g., gamma-ray bursts ) are particularly useful. 
In |2] we discuss scattering by a distribution of pebble 
sizes. We summarize our results in <J5] 

2. BASIC IDEA 

Here we show that small angle scattering of optical 
photons by large dust grains - "brilliant pebbles" - pro- 
vides a direct method of determining the size distribution 
of very large dust grains, and can also provide informa- 
tion on the spatial distribution of the dust. 

Particles with radii a large compared to the wavelength 
A have total extinction cross sections equal to twice the 
geo metric area (the so-called " extinction paradox" - see, 
e.g.. lBohren fc Huffman! Il983l ) : 50% of the extinction is 
contributed by diffraction of light, with much of the en- 
ergy concentrated in a forward scattering lobe. The char- 
acteristic angular radius for the scattering is 



A 



= 29'' 



A 



0.44/xm 



1 mm 



(1) 



For small scattering angles 9 S , the differential scattering 
cross section dC sca /d£l « constant, and for large angles 
dC sca ,/dfl oc 0~ 3 . It can be approximated by 



dC sl 



dfl 



1 + (1.80./0 O ) S 



(2) 



Figure [T] shows exact results for dC S c a /df l cal- 
culated for spher es using Mie theory (|Mid 119081 : 

iBohren fc Huffman! Il98l . Results are shown for radii 
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Fig. 1. — Differential scattering cross section vs. scattering an- 
gle for photons of wavelength A = 0.55 fim and spheres of radii 
a = 100 fim and 1 mm. The central lobe dominates the scatter- 
ing power in both cases. The results are insensitive to composi- 
tion: the figure actually shows results for both amorphous silicate 
spheres and graphite spheres, but the curves fall on top of one an- 
other. Also shown is cq. J2), which is seen to provide a very good 
approximation to the overall distribution of scattered power. 



a = 100 /im and 1 mm. Eq. ^ provides a good ap- 
proximation to the scattering for single spheres if one 
smoothes over the maxima and minima in the diffraction 
pattern. 

Assume that "pebbles" of radius a p account for a frac- 
tion T p of the integrated local mass density of dust pd 
where 
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da 



a,p(dn/da) p 



.{dn/ da) 



cha 
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Here, po is the density within a grain, and a c har ~ 0.15/im 
is the characteristic dust grain radius that is responsible 
for most of the dust mass and the bulk of the visual 
extinction in the Galaxy's interstellar medium. Along 
any given sight line, the ratio of visual optical depth of 
pebbles (with radius a p , and number ~ a p (dn/da) p ) to 
that of the entire grain population is given by 



At v (a p 



^char 



T p < 1. 



(5) 



Therefore such pebbles, if present, contribute negligibly 
to interstellar extinction. 

Now imagine a point source, that could be variable, as 
depicted in Figure El with an apparent visual magnitude 
my, located at a distance Do from an observer. Fur- 
thermore, assume that at distance D p from the observer 
there is a thin intervening dusty screen of large angular 
extent contributing a visual extinction Ay. 3 Scattering 
by the pebbles will produce a halo around the source 

3 GRB 050724's time-dependent X-ray halo shows that in- 
terstellar dust towar d this source is distributed in thin sheets 
HVauehan et al.ll2006ft . 



variable source 




Dn = D p /x p 

Fig. 2. — Geometry of optical light scattering by pebbles embed- 
ded in a dusty scattering screen. Examples of the variable source 
are novae, variable stars and GRB optical afterglows. 



with halo angle 9h ~ (1 — x p )6 s , where x p = D p /Dq; the 
characteristic halo angle is (1 — x p )8q. The time delay 
Ats due to small angle scattering is 



1 



D P e 2 h 



1 — x p 2c 
1.2 x 10 4 s 

1 — X T i 



(6) 
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1 kpc 



100' 



With the help of eq. ([5]), the photon flux scattered by 
pebbles of size a p is 



F s (t + At s )~F {t) 
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Note that in the above expression, Fo(t) is the "direct" 
photon power intercepted by the detector and Fg(t + 
Ats) is the photon power that underwent small-angle 
scattering off of pebbles, integrated over a scattering halo 
with radius 6% w (1 — x p )6 s with surface brightness 



(i\,s(t + At s ) 
mag arcsec -2 
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where m\(t) is the apparent magnitude of the point 



source. The solid angle AJ7/, 



of 



the scattering halo depends upon the size of the pebbles 
responsible for the small-angle scattering. With the help 
of cq. ([7]), we note that 
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Somewhat counter-intuitively, for fixed mass fraction T v 
the halo surface brightness increases with increasing peb- 
ble size a p . Even though the total power in scattered 
light goes as Fs oc a" 1 , the intensity increases oc a p be- 
cause the solid angle that the halo subtends decreases as 
a/. 

3. ASTROPHYSICAL APPLICATIONS AND FEASIBILITY 
OF DETECTION 

In the extreme event that the mass fraction of pebbles 
is large such that T p ~ 0(1), the scattered halo pro- 
duced by brilliant pebbles is still quite faint. By taking 
an extreme value of T p ~ 1 in cq. 0, we see that the 
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total power scattered is small such that F$(t + Ats) ~ 
10 -3 F (t). 

From Figure [1] we have = (1 — x p )6 s . Then, for dust 
in a single screen (i.e., single time delay for given halo 
angle 0) , the surface brightness of the halo core is 



Ma,, 



( (t + At s 



mag arcsec - ' 1 
m\(t) + 19.9 -2.5 log 



TyJ-p a p I /tmy 
(1 — Xp) 2 mm V A / 



(10) 



Of course, the above expression is independent of time 
for a source that does not vary. 

3.1. Strategy and Feasibility for Detection 

How does one detect a faint - possibly time variable - 
halo? 

For a p ~ 1 mm pebbl es, the scattering halo i s tens 
of arc seconds in radius. iDalcanton fc Bernstein! (|2000l ) 
obtained accurate surface photometry in the B- band 
where the sky is the most dim - of features as faint as 30 
mag arcsec -2 on scales larger than 10", using exposure 
times of 45 minutes on a 2.5m telescope. At a dark site, 
the limiting B surface brightness (for S/N=l) for a ~10" 
region will be 



:29.5 + 2.51og 10 (£>/2.5m) + 



(11) 



mag arcsec" 

1.25 log 10 (T/45 min) 

Thus A/ifl a t w 31 mag arcsec -2 could be achieved with a 
1 hr exposure on an 8 m telescope. The limiting factor 
will not be sky subtraction, but rather subtraction of the 
telescope PSF. 

The PSF of an optical telescope is typically charac- 
terized by a Gaussian "core" produced by atmospheric 
turbulence, a -3 halo resulting from diffraction by the 
large-scale telescope structure (e.g., aperture) at inter- 
mediate angles, followed by a "aureole" , varying roughly 
as -2 , due to a spectrum of small-scale i mperfections 
(e.g., micror i pples and dust o n the mirrors ) (lKing|ll971l : 
lRacine)[l99l lBernstein|[2007l ). 

The PSF of the Du pont 2.5m telesco pe at Las Cam- 
panas was measured bv lBernsteinl (|2007l ) in 2000 Septem- 
ber (~2 months after mirror realuminization) and can be 
approximately described at large angular radii by the fol- 
lowing broken power law 

7.5 log 10 (0/40") 6"<0<4O" (12) 



Mpsf ~ m 
« m 
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5.Olog lo (0/4O") 4O"<0<2OO" (.13) 

where /Lt ps f is in mag arcsec -2 . From here on, we use this 
exemplary PSF in our estimates. 

To estimate exposure times, assume that uncertainties 
in the sky brightness on 10" scales from an exposure time 
T sky is given by 

1/2 



AL 



sky 



= C 



1 



(14) 



sky \ ^sky-^sky , 

where C is a constant. Then, on 10" scales, for an ex- 
posure time T ps f, the uncertainty in the psf intensity at 
some angle will be given by 
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(15) 



Setting A/p S f = I p ^J(S/N), we can estimate the re- 
quired exposure time: 



T pst = T sky (S/N)< 



AL 



sky 



-'psf ~t~ -'sky 



(16) 



-^p.h. / \ -'sky 
~7sky(5'/^) 2 10^ ' 8( ' Alp " h " -AAtsky ' )-0 ' 4( ' Mpst-AlE ' ky ' l '(17) 

where we assume J ps f » 7 s k y . IDalcanton fc Bernstein! 
(|2000L 12002ft obtained A/i sky = 29.5 mag arcsec -2 for 
^ s ky = 22.2 mag arcsec -2 with T s k y = 45 min on the 
Dupont 2.5 m. Thus, if D is the telescope aperture, 

T psf = 0.75hr (^^j (S/N) 2 10°- 8 ^ -° A ^- 14 ' 72 

(18) 

is the exposure time necessary to detect the pebble halo 
with signal-to-noisc ratio S/N, assuming prior knowledge 
of the telescope psf. 

We now assess the detcctability of pebble halos. We set 
a p = 1 mm, Tp = 0.1, A = 0.44 /im, x p = 0.5, and ry = 
1 in following estimates. The characteristic scattering 
angle 0o = 29", resulting in a scattered halo with nearly 
constant surface brightness out to 6h = 14.5", with (from 
eq.QI)} 

Mp.h. = (m + 19.11) mag arcsec -2 (19) 



For the psf given by eq. I|12I13|) . at = 14.5", a steady 
source has 



%>sf = (m + 13.69) mag arcsec 
3.2. Statistical Detection 



(20) 



The column density of large dust grains may be cor- 
related with the column density of the sub- /im grains 
that produce interstellar reddening. It follows that a 
statistical search for brilliant pebble halos using large- 
scale imaging, such as the Sloan Digital Sky Survey, 
may be fruitful. Rather than fit for only an "intrin- 
sic" PSF Po(9,4>) that is the same for every star in 
a field (or at most dependent only on position within 
the field of view), one could determine whether the ac- 
tual images (excluding regions with small-scale nebular 
emission) can be better reproduced using a stellar PSF 
P Q {8,4>)+E(B - V)Pi(6), where 
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l \yh 
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E(B - V) 
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(1 - x p ) 2 V dft 



(21) 



is the dust contribution to the PSF, per unit E(B — V), 
where f(x)dx is the fraction of the reddening contributed 
by dust with x p in [x,x + dx] . For a first search, it would 
be adequate to assume the dust to be halfway to the star 
(i.e., f(x) — -> S(x — 0.5)). Using a large number of stars 
to determine Pi (0) may allow detection of the "brilliant 
pebble" phenomenon at levels that might be impossible 
for single stars. 

3.3. Bright Stars 

The pebble halo can be detected with reasonable ex- 
posure times on individual bright stars. For a star with 
Ay ~ 1 mag and tub < 8.5 mag (so that /i ps f < /Usky) eq. 
(TTSl) becomes 



-j- ) « 0.75hr(S/iV) 2 



2.5 m 
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10 



0.4m B -4.91 



(22) 
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Thus an exposure time of only 1.4 minutes would suf- 
fice for detection (with S/N = 5) of the pebble halo 
around a star with tub = 5 mag. The challenge will 
be to know the telescope psf well enough to accurately 
subtract it. This can be accomplished by measuring 
the psf using stars with minimal reddening. Because 
the uncertainty in the psf should be small compared to 
Ip.h./Ipsf ~ 0.007(Ay/mag), determination of the psf 
will require care. 



3.4. GRB Optical Afterglows and Other Bright 
Transients 

If the source is variable, the ability to detect a brilliant 
pebble halo of scattered light may increase dramatically. 
If the PSF of the telescope and atmosphere is stable and 
known, then point sources whose angular positions are 
fixed can be subtracted from each successive exposure. 
In doing so, the pebble halo, which varies in both time 
and space, may be isolated from the persistent stellar 
and interstellar background light. 

Gamma-Ray bursts (GRBs) occur at a rate ~ once per 
day in the Universe. Some GRB afterglows have opti- 
cal luminosities that can compete with those of quasars. 
For example GRB 990123 reached a peak optical flu x 
slightly brighter than my = 9 mag (jAkerlof et al.lfl999h . 
and GRB 0803 19B had my < 6 mag for ~ 40 sec 
(Bl oom etaLl f2008h . There may also be other ver 
bright optical transients: e.g.. lShamir fc Nemironl (|200 
observed what appeared to be a 5th magnitude flash 
that lasted ~ 10 min. In the event that a burst with 
rriB = 5 mag were to occur at low Galactic latitudes, the 
rapid fading of the point source provides the ideal circum- 
stance for detection of a pebble halo, because observa- 
tion of the time-delayed scattered-light halo can be done 
without interference from the telescope psf. If an opti- 
cal transient shone with mg = 5 for a few minutes and 
then faded, the time-delayed scattered light "halo" (with 
\i = 24.1 magarcsec" 2 if dust with Ay = 1 is present in 
a thin sheet) could be detected in ^2 min of integration 
on an 8 m class telescope. 

3.5. Variable Stars 

There are several thousand variable stars in the M ilky 
Way brighter than my = 10 (jPaczvhski et al.ll2006l ). In 
fact, from the ^45^45 catalogue, there are over 400 eclips- 
ing binaries with periods shorter than 1 day and V < 10. 
For ~ 100^m pebbles the delay timescale is in terms of 
hours, rather than tens of minutes. It therefore follows 
that the pebble halo can be detected by image subtrac- 
tion if the binary is significantly reddened and distant, if 
~ 100/im grains are sufficiently numerous. 

4. SMALL ANGLE SCATTERING FROM A SIZE 
DISTRIBUTION OF PEBBLES 

So far, we have considered small angle scattering due 
to the presence of pebbles of a single characteristic size. 
Now, we briefly consider a distribution of sizes. For sim- 
plicity, we assume that the pebble size distribution re- 
sembles a power law 



dN 
da 



Afa 



< a < a r 



(23) 



where N(a) is the column density of pebbles with radii 
a p < a. For p = 4 (i.e., equal mass per logarithmic 



interval in a p ) the normalization constant is given by 

T V Ochar 



Af 



(24) 



2n ln(a max /a n 

where the optical depth Ty is assumed to be provided by 
the dominant dust population with radii a c h ar ~ 0.15/xm, 
with extinction cross sections C ex t ^ 7Ta \±_ a , r - 

For scattering by material in one scattering screen, we 
define the differential scattering optical depth such that, 
at halo angle 9) ll the scattered intensity is 

F (t) ( dr., 



I{0 h ,t + At s ) 



where 
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(25) 



(26) 



9 s =e h /{\-x v ) 

With the help of the approximation {2J for dC SCSi /dn, the 
differential scattering optical depth is 

dH A 2 J atnin l + (a/a c ) 3 1 ' 

(1_^)A 

1.8n0 h V ' 

We can distinguish three regimes: the core, with Oh < 
6i = (1 — x p )A/1.87ra max ; the intermediate halo, with 
#i < Oh < #2 = (1 — Xp)X/1.8ira m i n ; and the outer halo, 
with 0^ > 6 2 ■ For 2 < p < 5 we have 
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(29) 
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(5 - p){p - 2)(1.8tt) 5 ~p \p~ 3 5 h ~P 

for 6 l < 9 h < e» 2 (30) 
tt 2 (1-xJ 3 A 



for h > 2 (31) 

In Figure [3J we perform the integral in eq. (|27[) over 
a p = 4 pebble size distribution for A = 0.4405 fiui and 
0.802 fjja. As an illustrative example, consider the case 
in Figure [3] where the large pebble size cutoff a max = 
1.0 cm for A = 0.4405 /im. The angle 0\ marking the 
transition from the the core, where dr s /dQ ~ const, to 
the intermediate halo, where dr s /dQ oc h ~ v , indicates 
the value of the upper cutoff in pebble size, a max ~ (1 — 
a; p )A/1.87T0i. 

The power-law index p for the size distribution can 
be determined from both the angular dependence 0^ -5 
and the color (dT s /dfl oc A 3 ~ p ) of the intermediate halo 
(6*i <C 0h <C #2). The large angle behavior dr sca /df2 oc 
0~ 3 for > 2 indicates that at these angles, all of the 
pebble size range [a m i n , a max ] is in the large scattering 
angle regime. 

Figure [3] shows how dr sca /df2 varies with photon wave- 
length. For shorter wavelengths, the halo makes the tran- 
sition from the small angle to intermediate angle regime 
at smaller angles. The scattered halo will be very blue 
in the core: for p — 4 the core will have a color excess 
B - I = -2.51og 10 (0.8020/.4405) 2 = -1.30 mag rela- 
tive to the unscattered source, while the intermediate 
zone beyond the core will have B — I = —0.65 mag. For 
> 02 the halo will be red, with B — I = +0.65 mag, but 
the combination of long delay time Ats and low surface 
brightness may render this outer halo undetectable. 
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0.1 1 10 100 1000 10" 



0s =6y(l-x p )[arcsec] 

Fig. 3. — Pebble size-averaged differential optical depth 
(^Tsca / dQ) I (ry T p ) vs. scattering angle 8 a for photons of wave- 
length A = 0.4405 (im (blue) and A = 0.802 (im (red). Each curve 
represents dC BC a/dQ, integrated over an a -4 pebble size distribu- 
tion, but with varying ranges in pebble size. The lower limit for 
pebble sizes a m i n = 10 fim for all of the curves. The maximum 
pebble size a max = 1 cm, 1 mm and 10 2 /im for the top, middle 
and bottom curves, respectively, for each of the two photon wave- 
lengths. At small scattering angles, dT BCB ,/dCl oc a constant, while 
dT S ca,/dQ oc 6s~ 3 at large scattering angles, recovering the asymp- 
totic behavior of eq. (2}. At intermediate angles, dT sca /df2 oc . 
This 'intermediate" angle range is more prominent for broad dis- 
tributions of pebble sizes. The dashed lines at the bottom repre- 
sent the limit where the pebbles are replaced purely by "normal" 
a ~ 0.15 (im grains for each photon band. As long as T v > 10 — 2 for 
the pebble-size distributions utilized in this figure, then the vari- 
able scattered light should be dominated by pebbles rather than 
normal dust. 



at a distance ~ 1 kpc. We refer to this phenomena as 
"brilliant pebbles." 

Detection of the brilliant pebble halo will be technically 
challenging. Even for favorable assumptions regarding 
the pebble size distribution, an accurate determination 
of the telescope PSF is required in order to detect such a 
low surface brightness signal. Nevertheless, we show that 
a pebble halo could be detected using existing facilities 
if the mass in ~1 mm pebbles is more than a few percent 
of the total dust mass, as has been suggested by observa- 
tions of high- velocity dust grains in interplanetary space 
and hyperbolic micrometeors (see § 1). 

Detection is best accomplished using a wcll-bafflcd tele- 
scope with a freshly aluminized mirror, to minimize scat- 
tering by small-scale imperfections such as dust on the 
mirror. 

A brief but bright optical burst would be ideal, as the 
direct light can fade while the time-delayed pebble halo 
persists, allowing the telescope PSF and the pebble halo 
to be separated in the time domain. Also, the pebble 
halo from a bright short period eclipsing eclipsing binary 
that is significantly reddened can be subtracted from the 
PSF in the time domain as well. 

If a "brilliant pebble" halo is detected, the angular 
variation of the halo intensity will indicate the pebble 
size range and the form of the particle size distribution. 
The halo will be very blue in the core, and moderately 
blue in the intermediate halo; the color of the interme- 
diate halo provides an independent measure of the size 
distribution. If the source is suitably variable, the dis- 
tance to the pebbles can be determined from the time 
delay of the halo relative to the point source. 

Upper limits on brilliant pebble halos can provide valu- 
able constraints on the size distribution of solid particles 
in the interstellar medium. 



5. SUMMARY 

We propose a method for the detection of extremely 
large (~ 0.01 — 1 cm) dust grains. When the ratio of pho- 
ton wavelength to grain size A/a <C 1, scattering of light 
is highly peaked in the forward direction, with a char- 
acteristic scattering angle r?o ~ A/Va ~ 30(1 mm/a p )". 
For a variable source, the halo will lag with a character- 
istic time-scale ~ 10 min for scattering by 1 mm pebbles 
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